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Abstract 

We explore the phenomenological consequences of general late-time modifications 
of gravity in the quasi-static approximation, in the case where cold dark matter is 
non-minimally coupled to the gravitational sector. Assuming spectroscopic and pho¬ 
tometric surveys with configuration parameters similar to those of the Euclid mission, 
we derive constraints on our effective description from three observables: the galaxy 
power spectrum in redshift space, tomographic weak-lensing shear power spectrum and 
the correlation spectrum between the integrated Sachs-Wolfe effect and the galaxy dis¬ 
tribution. In particular, with ACDM as fiducial model and a specific choice for the 
time dependence of our effective functions, we perform a Fisher matrix analysis and 
find that the unmarginalized 68% CL errors on the parameters describing the modi¬ 
fications of gravity are of order a ~ 10“^-10“^. We also consider two other fiducial 
models. A nonminimal coupling of CDM enhances the effects of modified gravity and 
reduces the above statistical errors accordingly. In all cases, we find that the parame¬ 
ters are highly degenerate, which prevents the inversion of the Fisher matrices. Some 
of these degeneracies can be broken by combining all three observational probes. 
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1 Introduction 


The recent measurements of the cosmic microwave background (CMB) anisotropies, per¬ 
formed by the WMAP and Planck satellites, have significantly improved our knowledge 
on the content of the universe and on the initial conditions of cosmological perturbations. 
A similar progress is expected from the next generation of galaxy surveys concerning the 
properties of dark energy or, possibly, modifications of general relativity on cosmological 
scales. Indeed, even if the CMB is useful to constrain dark energy through the integrated 
Sachs-Wolfe (ISW) effect and gravitational lensing, these effects are ultimately related to 
the impact of dark energy on the late-time evolution of structures. Probing directly these 
large scale structures is thus thought to be the most promising source of information on 
the origin of the current acceleration. 

Since no compelling model of dark energy has emerged from theoretical investigations, 
it is appropriate to resort to a description that encodes a wide range of physical effects 
with a limited number of theoretically motivated parameters, in order to compare devia¬ 
tions from the standard ACDM scenario with cosmological observations on linear scales. 
For single-held dark energy models in the presence of universally coupled matter helds, 
this research program has been initiated by the effective theory of dark energy recently 
proposed in Refs. [1-3], inspired by the so-called effective held theory of inhation [4,5] 
and of minimally coupled dark energy [6]. Another model-independent framework that 
has been developed with the same motivations is the Parameterized Post-Friedmann ap¬ 
proach [7,8] . In the effective theory of dark energy, the quadratic action describing linear 
perturbations of single-held models belonging to Horndeski theories is characterized by 
four free functions of time [3,9-11], while a hfth function must be introduced to describe 
theories beyond Horndeski [12,13]. The power and efficiency of this formalism has just 
started to be exploited. For instance, it has been applied to explore and forecast the 
phenomenology of dark energy and modihed gravity in [14-17] (see also [18,19] for some 
nonlinear aspects). 

Recently, in Ref. [20] , we extended this unifying treatment to allow for distinct confor¬ 
mal and disformal couplings of matter species to the gravitational sector.^ We focused on 
Horndeski-like models, i.e. those whose quadratic action has the same structure as Horn¬ 
deski theories,^ although the full action can be different. This is a rather natural extension 
given that a modihcation of the gravitational sector can often be interpreted as a direct 
coupling of matter to a hfth force exchanged by the scalar, in the frame where the scalar 
and the gravitational huctuations are demixed—the so-called Einstein frame. Together 
with the four functions describing the gravitational quadratic action, each matter species 
is now characterized by two new functions parametrizing their conformal and disformal 
couplings to the gravitational metric. However, as reviewed in Sec. 2, the structure of 
the full action remains invariant under conformal and disformal transformations of the 
gravitational metric itself. Taking into account this freedom, which allows for instance to 

treatment of single-field dark energy coupled to CDM in the context of the Parameterized Post- 
Friedmann framework can be found in [21]. 

^Note that although Horndeski theories are generically unstable under quantum corrections [22], an 
example of a radiatively stable subclass of Horndeski theories where all the operators of action (2.3) can 
be relevant has been proposed in [23], based on weakly broken galileon invariance, and applied to inflation 
in [24]. 
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choose a frame where one of the species is minimally coupled, one eventually finds that the 
whole system depends on a total of 2{Ns + 1) independent functions of time, where Ns is 
the number of matter species. In this context, the conditions for stability (i.e. the absence 
of ghostlike and gradient instabilities) can be generalized to any frame (see Sec. 2). 

In this article we go one step further and explore the constraining power of future large 
scale structure surveys on the deviations from the standard ACDM scenario, expressed in 
terms of the parameters of the effective theory of dark energy proposed in [20] . Specifically, 
we will consider a simple scenario where the gravitational sector is described by Horndeski¬ 
like models while, in the matter sector, cold dark matter (CDM) is nonminimally coupled 
to gravity. This extends to a much broader spectrum of gravitational theories previous 
studies of coupled dark energy, with conformal [25,26] (see also [27] and references therein) 
and disformal (see e.g. [28-37]) couplings. 

The equations of motion for the linear perturbations in the presence of modified gravity 
and nonmininally coupled CDM, derived in [20], are reviewed in Sec. 3, where we assume 
the quasi-static approximation. As shown in [38] , this approximation should be reliable for 
surveys such as Euclid as long as the sound speed exceeds 10% of the speed of light, i.e. Cg > 
0.1. In particular, we will consider the extreme quasi-static limit, i.e. the limit k —>• oo, of 
the dynamics. In such a regime the linear growth of matter (both for baryons and CDM) 
remains scale-independent as in ACDM. Modifications of gravity and the nonminimal 
coupling to CDM are encoded in the time dependence of the gravitational couplings in the 
“Poisson” equations for the metric potentials, which are different for baryons and CDM. 
As explained in Sec. 3, this time dependence modifies the growth rate of structures and 
the lensing potential, which in turn affect, respectively, the redshift-space distortions and 
the weak-lensing cosmic shear. 

In Sec. 4 we introduce the details of our parametrization, in particular concerning 
the time dependence of the parameters characterizing the modifications of gravity. We 
consider three fiducial models: a minimal ACDM model, a braiding model and a model 
with an active nonminimal coupling of CDM. In Sec. 5 we perform a Fisher matrix analysis 
based on future photometric and spectroscopic data with configuration parameters close 
to those of the Euclid mission [27,39] as an example. We focus on the two-point statistics 
and consider the galaxy power spectrum in redshift space for the spectroscopic data, the 
projected weak-lensing shear power spectrum for the photometric data as well as the 
correlation between the ISW effect in the CMB temperature and the photometric galaxy 
distribution. The derived constraints are discussed in Sec. 6, together with the involved 
degeneracies. It should be mentioned that other approaches have been developed to study 
in a general and model-independent way the impact of modified gravity on cosmological 
observables, together with the involved degeneracies, e.g. on the growth rate of fluctuations 
[40] (see also [41,42]) or on the weak lensing [43]. 

In Sec. 7 we summarize our results and draw conclusions. Details on the parametriza¬ 
tion and the choice of background cosmological parameters are given in the App. A, while 
in App. B, we discuss the frame dependence of the evolution equations of matter. 
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2 Model and main equations 


In this section, we introduce our general formalism and then focus on the specihc model 
at the core of the present work. The first subsection, which is mainly a review of our 
recent paper [20] and previous works, can be skipped by the reader mostly interested in 
our phenomenological model and forecasts for the parameter constraints. The model that 
we are specifically studying in the rest of this paper is described in the second subsection. 

2.1 Effective description of the gravitational and matter sectors 

We start by summarizing the effective approach of dark energy introduced and developed 
in Refs. [1,3,20] (see e.g. [11,44] for reviews). The gravitational sector is assumed to consist 
of a four-dimensional metric and of a scalar field cf). In order to treat simultaneously 
a wide range of models, it is very convenient to “hide” the scalar field in the metric, 
by choosing the constant-time hypersurfaces to coincide with the uniform scalar field 
hypersurfaces. In this gauge, referred to as unitary gauge, the metric can be written in 
the ADM form [45], 

^ + hij (dx* + N^dt) {dx^ + NHt) , (2.1) 

where N is the lapse function, N'‘ the shift vector and hij the three-dimensional spatial 
metric. 

In unitary gauge, a generic gravitational action can be written in terms of geometric 
quantities that are invariant under spatial diffeomorphisms, namely in terms of the lapse 
N, the 3d Ricci tensor Rij of the constant time hypersurfaces, as well as their extrinsic 
curvature ATjj, with components 

~ ~ “ ^j^i) ) ( 2 - 2 ) 

where a dot stands for a time derivative with respect to t, and Di denotes the covariant 
derivative associated with the spatial metric hij (spatial indices are lowered or raised via 
the metric hij). 

The generalized Friedmann equations are then obtained by varying the specialization of 
the action to a homogeneous FLRW (Friedmann-Lemaitre-Robertson-Walker) spacetime, 
endowed with the metric = —dt^ + a?{t)dx^ . The dynamics of the linear perturbations 
is governed by the quadratic action, obtained by a perturbative expansion of the original 
action. 

In this paper, we will consider a very large class of models, which includes all Horndeski 
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theories, for which the quadratic action can be written in the form [3,9-11]^ 



,3 , 

d xdt ^ 


6K}6Kl 


6 K^ + R5N + (1 + ot) <52 (VhR/a^^ 


+ aKH^6N^+ 4aBH6K6N , 


(2.3) 


where M, ot, Q^b and ax are four time-dependent functions and 82 denotes the second 
order term in a perturbative expansion. H = d/a is the Hubble parameter. We have 
not included irrelevant terms that vanish when adding the matter action and imposing 
the background equations of motion. Note that (2.3) does not include the models beyond 
Horndeski [12] for which the coefficient of the term R 6 N differs from 1, the difference 
defining a new parameter an [Hj- 

General relativity corresponds to the particular case where ot = ob = OiK = 0 and 
M = Mpi. In general, the above quadratic action contains not only two tensor modes, 
as in general relativity, but a scalar mode as well. The coefficient in front of the tensor 
kinetic term is and, by analogy with general relativity, M can be identified with an 
effective Planck mass. If M depends on time, it is convenient to introduce the related 
parameter 


CTm = 


1 dlnM^ 
H dt 


(2.4) 


The parameter appears in the gradient term of the tensor modes and is thus directly 
related to the tensor propagation speed, namely 


= 1 + aT- 


(2.5) 


The stability of the tensor modes is ensured by requiring > 0 (absence of ghosts) and 
aT > — 1 (absence of gradient instabilities).^ 

Keeping in mind that the lapse perturbation is analogous, in the ADM language, to the 
time derivative of the scalar perturbation, one observes that the parameter ax is related 
to the coefficient of the kinetic scalar term. It is thus present for simple quintessence 
models. Finally, the coefficient ob characterizes the mixing between the scalar and tensor 
kinetic terms, sometimes called “braiding”. In contrast with the tensor modes, the full 

^Together with the operator a^iSNR, this is the most general quadratic action for linear perturbations 
about a homogeneous and isotropic spacetime that does not induce higher derivatives in the equation 
of motion of the linearly propagating scalar degree of freedom. In consistent effective theories, higher 
time derivatives are not forbidden but are suppressed by positive powers of the ratio between the energy 
and the cutoff scale (see e.g. [46,47]). Thus, at energies much smaller than the cutoff their effect can be 
neglected without loss of generality. Higher spatial derivatives are not necessarily suppressed and may 
dominate the dispersion relation, such as in the Ghost Condensate theory [48]. In this case, higher spatial 
gradients become relevant, and can easily be included in our formalism, but begin operating at very short 
distances [6,49], typically shorter than the cosmological ones. 

'^As shown in ]20,50] and reviewed below, the propagation speed for gravitons can be set to unity by a 
convenient disformal transformation (only ratios between sound speeds are invariant and thus meaningful 
physical quantities). It is thus not a priori pathological to have Cy > 1 in a generic frame and we will 
not impose any upper bound on ct as a condition for the viability of the theory. A propagation speed 
for gravitons smaller than that of the other particles is instead very tightly constrained at high energy by 
cosmic rays observations [51]. We have not taken this bound into account in our analysis, since it concerns 
the speed of gravitational waves at wavelengths much shorter than the cosmological ones. 
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dynamics of the scalar mode depends on the matter action as well, and the discussion on 
the scalar stability conditions thus needs to be postponed until after the introduction of 
the matter action below. 

The remarkably simple form of the quadratic action (2.3) holds only in the unitary 
gauge. However, it is straightforward to derive the quadratic action in an arbitrary gauge, 
by simply performing a time reparametrization of the form 

t —)>(/> = t + Tr{t, x) , (2.6) 

where the unitary time becomes a four-dimensional scalar field. The scalar degree of 
freedom of the gravitational sector thus reappears explicitly in the form of the scalar 
perturbation vr. 

A matter species can be either minimally or nonminimally coupled to the gravitational 
metric g^u- In the latter case, it is often assumed that matter is minimally coupled to 
some effective metric which depends on and on the scalar field (f). We will adopt 
this type of nonminimal coupling in the following and consider a matter action of the form 


*5111 — *5m [V’ni) 5/11* 


(2.7) 


with 

gfiu = C{(j))g^u + . ( 2 . 8 ) 


The initial gravitational metric gf^i, being somewhat arbitrary in general, one has 
the freedom to choose the metric as the new gravitational metric. Remarkably, the 
quadratic action (2.3) remains of the same form [20,52], ^ with new parameters defined as 




m2 

C*\/l + ctd 


(2.9) 


and® 


OK 


OiB 


au 

Ctq' 


Ok + 12aB[aD + (1 + aD)ac] ~ 6 [q;d + (1 + aD)ac]^ + SHmOD 

(1 + ac)2(l + aD)2 

1 + CkB ^ 

(1 + Q!c)( 1 + cto) 

ttM — 2ac «D 

1 + ac 277(1 + aD)(l + ac) ’ 

(1 + aT)(l + ao) — 1 , 


( 2 . 10 ) 


where 


ac = 


C 

2HC ’ 


ao = 


D 

C-D' 


( 2 . 11 ) 


Given a single species of matter, one can thus always work in the frame where this species 
is minimally coupled. If there are several matter species, this is possible only in the case 


®In the presence of the operator proportional to qh [3, 11] describing linear perturbations in the 
theories beyond Horndeski proposed in [12,13], the structure of the Lagrangian remains invariant under 
the transformation (2.8) even if the disformal function D depends on (9</>)^ as well [12] (see also [53] for a 
recent study). 

®Here we correct a typo in the expression for Sk in eq. (2.45) of the arXiv version of Ref. [20]. 
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of universal coupling, i.e. if all species are coupled to gravity via the same effective metric. 
By contrast, for species with different couplings, one cannot find a frame where all of them 
are minimally coupled. It remains however possible to choose a frame where one of the 
species is minimally coupled, even if the others are not.' 

The sum of the gravitational and matter actions at quadratic order yields the dynamics 
of the scalar mode, as mentioned earlier. As shown in [20], the kinetic term of the scalar 
mode is proportional to the combination 


Oi — ckk 


+ GfUg + 3 'y ^ ctd,/ 


where 


0/ = 


PI 




( 2 . 12 ) 

(2.13) 


while its propagation speed is given by 


Cg — — — { (I+ob) 


-jp—OiM+OiT+a-Qi'^+OiT:) 


ctR 3 
H—--\— 

H 2 


|^l+(l+aD,7')?i’7' • (2.14) 


The stability conditions for the scalar mode. 


a > 0 , > 0 , 


(2.15) 


involve all the modified gravity parameters, as well as the matter disformal couplings. 


2.2 Baryon-CDM model 


In our model, the coupling of CDM to the gravitational sector is different from that of the 
other species (baryons, photons and neutrinos). In the following, for simplicity, we choose 
to work in the frame where the other species are minimally coupled and assume that the 
original metric corresponds to this frame (if not, one just needs to apply the above 
metric transformation). We then assume that the coupling of CDM to gravity and dark 
energy is characterized by an effective metric of the form 

, (2.16) 

from which one can define, in analogy with (2.11), the conformal and disformal parameters 


ac,c 


Cc 

2HCc ’ 


ao.c 


Dc 

Cc-Dc' 


(2.17) 


We ignore the photon and neutrino cosmological fluids, as we are interested in late-time 
cosmology where their effects are negligible. 

The equations of motion for the matter species follow from the conservation, or non¬ 
conservation, of their respective energy-momentum tensor. Since baryons are minimally 
coupled, their energy-momentum tensor is conserved as usual, i.e. 


V^r(fe)^, = 0. (2.18) 

^The situation simplifies during inflation, when the couplings to matter can be ignored. In this case, 
without loss of generality one can always go to a frame where ayi = dx = 0, corresponding to the standard 
time-independent Planck mass and unity speed of propagation for gravitons. In this frame one then recovers 
the standard inflationary predictions [50]. 









By contrast, the CDM energy-momentum tensor is not conserved, but instead satisfies the 
equation 


+ Qcdyfj) — 0 


(2.19) 


with 




where a prime denotes a derivative with respect to cf). Like the usual conservation equation, 
this equation can be derived by simply using the invariance of the matter action under 
arbitrary diffeomorphisms. 

The background evolution equations for the baryon and CDM fluids follow directly 
from (2.18) and (2.19). On a FLRW background, the definition of Qc, eq. (2.20), reduces 
to 


Qc 


Hpc 
1 + CKD,C 


|ac,c + aD,c 



dD c 1 

^ 2^/(1 +an,c) I ■ 


( 2 . 21 ) 


Substituting the above expression into eq. (2.19), one finds that the homogeneous fluid 
equations can be written in the form 


Pb + 3Hpb = 0 , 

Pc + 3H{1 - 7c)/0c = 0 , 


where the coupling parameter 7 c is given by® 


1 

7c = + 


dD,c 

6iL(l -h aD,c) ' 


( 2 . 22 ) 

(2.23) 


(2.24) 


Expressed in terms of the energy density fractions defined in (2.13), the evolution 
equations for the baryon and CDM energy densities, (2.22) and (2.23), become 


7 = -^^(3 + 2^ 

i2c - + 2j^ 


+ OiMj : 

— 37c -|- q;m 1 Dc . 


(2.25) 

(2.26) 


The presence of the coefficient om is due to the fact that the mass M, which appears in 
the definition (2.13), can be time-dependent. 

The evolution of the Hubble parameter is usually determined by the Friedmann equa¬ 
tions. In the present work where dark energy remains unspecified at the background level, 
one can alternatively assume some specific evolution H = H{t) and infer from it the dark 
energy background components. This means that the Friedmann equations, written in the 
form 

+ ^De) , H = [Pm + (1 + W^De)/5De] , Pm ^ Pb + Pc , 

(2.27) 

^Taking into account eq. (2.23) one finds that Qc = SHpcjc- 
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are treated as definitions of the energy density for dark energy, /Ode, and of its equation 
of state parameter, u^de, namely 


_2 H 

Pde = — pra , tfDE = t ^ n —' ’ (2.28) 

1 — i irn 

where 

n^ = nb + ^c- (2.29) 

Given some prescription for the time-dependent functions H = H{t), auit) and 7c(t), 
the evolution of Qi, and Qc can be determined in terms of their present values and 
Qc,o- This will be done explicitly in Sec. 4.1. 

3 Linear perturbations 

In this section, we present the equations governing the linear perturbations. For con¬ 
venience, we work in the Newtonian gauge, where the scalarly perturbed FLRW metric 
reads 

ds'^ = -(1 + 2 ^)dt^ + a‘^{t){l - 24 ')df 2 . ( 3 . 1 ) 

For each species, the continuity and Euler equations can be derived from, respectively, the 
time component and the space components of eqs. (2.18)-(2.19). As obtained in [20], they 
read in Fourier space 

h - -^Vb = , (3.2) 

Vb = -^, (3.3) 

^2 

Sc - -^Vc = 3(4' + jcHTTy + 2(1 -h aD,c)(ac,c - 37c)i4(4> - tt) - aD,c(^ - , (3.4) 

Vc + ^HjcVc = -4> - 3Fl7c7r . (3.5) 

These equations must be supplemented by the generalized Einstein equations and by 
the scalar fluctuation equation. We will not write them explicitly here but they can be 
found in [20]. 

3.1 Quasi-static approximation 

The evolution of perturbations well inside the horizon is most conveniently studied within 
the quasi-static approximation. This is justified for spatial scales that are smaller than 
the sound horizon of dark energy, or equivalently for wavenumbers k ^ aH/cs (see [38] 
for a detailed discussion and [54] for a recent analytical extension of this approximation). 
In this regime, one can neglect time derivatives with respect to space derivatives and the 
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continuity and Euler equations (3.2)-(3.5) for the baryon and CDM fluids simplify into 


^2 

- KVb = 0, (3.6) 

Vb = -^, (3.7) 

r.2 

4- 2^c = 0, (3.8) 

Vc + 3H-fcVc = -$ - 3i77c7r . (3.9) 


The equations for the gravitational potentials and T and for the scalar fluctuation 
TT also simplify and become constraint equations. The gravitational potentials satisfy two 
Poisson-like equations, given by [20] 

r.2 o 

— {(l + ckt + /?!) bObSb -|- [1 -I- ax + Wc^c} , (3.10) 

-2~ {(1 -|- /^b/?^) + [1 + /3b(/^5 + l^'y)] j (3-11) 

Qj Li 

where we have introduced the parameters a;/ = flz/rim, 


/I - ^ 

PB = -T77“B , 


1 /2 


\/2 \/2 

^5 ^ ^ „ „i/9 [“b(1 + ax) + ax - an] , 


(3.12) 


CsOL 


CgOL 


1/2 


as well as® 


/^7 — _ _ 1 /9 "Ic 


The scalar fluctuation also satisfies a Poisson-like equation, which reads 

o LJ-o + ^7 )^c' 5 c 

a"* \' 2 csa^/^ 


(3.13) 


(3.14) 


Combining eqs. (3.6)-(3.9) with eqs. (3.10)-(3.11) and (3.14) leads to a system of two 
second-order equations for the density contrasts. 


6 b + 2H6b — -H^Qm {(1 + + I6‘f)i^b6h + [1 + Q^t + /3g(/3g + Pry)] Wc<5c} j (3.15) 

+ (2 + 3jc)H6c = -H^^m {[1 + Q^T + /3^(/3g + ^b6b + [l + + P'y)'^] ^c6c\ ■ 

(3.16) 


Introducing the bias be (6b) between CDM (baryons) and the total matter density contrast 
(5m — 

6c = 6c 6m (6b = 6b 6m), (3.17) 

®The parameter P-y generalizes the parameter P defined for coupled quintessence in Sec. 5.3.4 of [15]. 
In this case, the relation between the two parameters is P-/ = —%/2/3. We thank Valeria Pettorino for a 
discussion on this issue. 
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the influence of modified gravity and nonminimal coupling onto the growth of perturba¬ 
tions enters through the combinations 


Tb — ot + + f3-fUJcbc) , Tc = ctT + (/5^ + l3'y){/3^ + P^^cbc) ■, (3.18) 

which vanish for standard gravity (the friction term 7 c on the left hand side of eq. (3.16) is 
essentially a background effect and does not affect directly the energy density perturbations 

bpbyC) ■ 

Modifications of gravity exchanged by vr are parametrized by (3^ and the nonminimal 
coupling of dark matter is parametrized by [20]. This separation of effects is not 
physical and depends on the choice of frame. Indeed, under a generic change of frame 
(2.8), one finds, using (2.9)-(2.10) as well as the relations 


«D,7 


an,I — od 
1 + aj) 


OiC,I 


ac,i - ac 
l + ac 


that these two parameters transform as 


+ /57*)(1 + , 

/37 = iP'y - / 37*)(1 + , 


(3.19) 


(3.20) 


where 


/^7* 


3\/2 


, do 
2F(l + aD) 


(3.21) 


See also App. B for a discussion on the frame dependence of eqs. (3.15) and (3.16) and of 
the combinations Tb^c- 

The modification of gravity associated with the parameter ctt does not depend on the 
exchange of vr, see eq. (3.14) and Refs. [20,55] (see also [56] for a recent discussion on 
local constraints of this effect), and does not mix with the other two effects under change 
of frame. We note that if > 0 (which corresponds to a speed of graviton fluctuations 
ct > 1) in the absence of nonminimal coupling, i.e. f3ry = 0, the combinations (3.18) are 
always positive, which tends to enhance the growth of structure. More generally, for a 
positive ot the combinations and Tc can be negative only if (3^ has the opposite sign 
of (3^. 

Since equations (3.15)-(3.16) are independent of the wavenumber k, one can factorize 
the time dependence from the k dependence of the initial conditions and write the solutions 
in the form 

6c{t, k) = Gc{t) 5cfi{k) , 5b{t, k) = Gb{t) Sbfi{k) , (3.22) 


where (5c,o and (5fe,o represent the initial density contrasts for CDM and baryons respectively, 
defined at some earlier time in the matter dominated era. The two functions of time Gc{t) 
and Gb{t) are the growth factors for CDM and baryons, respectively, assumed to be equal 
at the initial time, Gc(0) = Gfe(O) = 1. 

The continuity equation (3.8) then implies that the velocity potential Vc for CDM is 
given by 


Veit, k) 


a 

kP' 


Gc{t)5c,m 


a‘^H 

kP 


feit) 6cit, k) , 


(3.23) 


12 










where, in the second equality, we have introduced the CDM growth rate 


^ din Gc 
^ din a 


(3.24) 


Similarly, using the continuity equation (3.6), one finds that the velocity potential Vh for 
baryons is given by 


Vb = ^f{t)dh{t,k), (3.25) 

d In a 

3.2 Link with observations 

We now examine how the quantities introduced above can be probed by cosmological 
observations. 

A powerful cosmological probe for dark energy is weak leasing, which depends on the 
so-called scalar Weyl potential, i.e. the sum of the two gravitational potentials and 'k. 
Combining the Poisson-like equations (3.10) and (3.11), one gets the expression 

4> -|- 'k =- 2^2 '^b) + P'y^ckc)] <^m • (3.26) 

In analogy with the combinations (3.18), it is convenient to define 

Tl'lens ^ CUT + + Pb) iP^ + P^CJcbc) , (3.27) 

which vanishes when gravity is standard. 

Another way to probe dark energy is via the observation of galaxy clustering. In 
particular, redshift-space distortions are sensitive to the growth rate of fluctuations, which 
is affected by deviations from standard gravity. Here we extend previous studies and 
include also the effect of a nonminimal coupling of CDM. 

When observing galaxies, one must take into account the fact that what is directly 
measured is the redshift, and not the distance of the galaxy. In the parallel plane ap¬ 
proximation, the correspondence between the so-called redshift space and real space is 
described by the change of coordinates (see e.g. [57]) 

s = (3.28) 

an 

where s and x denote the spatial coordinates in redshift and real space respectively and 
Vg^z is the line-of-sight component of the galaxy’s peculiar velocity. At linear order, the 
invariance of the number of galaxies yields the expression for the number density in redshift 
space in terms of the number density in real space, 

bg,s = bg — • (3.29) 

On large scales, the galaxy peculiar velocity Vg can be related to the CDM and baryon 
fluid velocities, respectively Vb and Vc, by effectively treating galaxies as test particles 
(see e.g. [58]) of baryon and CDM mass fractions Xb = Mb/Mg and Xc = Mc/Mg (Mg = 
Mb + Me), respectively. By considering that the large-scale galaxy momentum coincides 
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with the sum of the baryon and CDM fluids momenta in the linear regime, the galaxy 
peculiar velocity is given as 

Vg = XcVc + XbVb , (3.30) 

where Vc = Vuc and Vb = Vvb are the linear velocities satisfying the Euler equations (3.7) 
and (3.9). Indeed, in the absence of screening the mass of the CDM component in the 
galaxy is not conserved and obeys 


Me = SHjcMc , 


(3.31) 


in agreement with the background evolution (since Me scales as PeO?)- Then, the combi¬ 
nation of the Euler equations yields 

— (MgVg) = — (MeVc) + ^ {MbVh) = Fg, (3.32) 

where 

Fg = —MgV'h -|- SHjeMe'VTT (3.33) 

is the neat force exerted on each galaxy. The last term is due to the fifth force on the 
CDM component. 

Using the expression (3.23) and (3.25) for the velocity potentials, one thus finds 


V 


g ~ 




{Xefe^c T Xbfb^b) ■ 


(3.34) 


Substituting the above expression into (3.29), and proceeding as in the standard calcula¬ 
tion, one hnally obtains, in Eourier space, 

k 

<5g,s — p {Xefe^c T Xb/b^b} j p — ~j^ > (3.35) 

or 


— ^g T t {Xefebe T Xbfbbb)bg, (3.36) 

"g 

after introducing the galaxy bias 6g, defined by 

Sg = bg5m- (3.37) 

The galaxy power spectrum in redshift space is thus given by 

PgAk) = Al + pAesfPraik) , (3.38) 

where we have introduced the effective growth rate of the galaxy distribution as 

feS = Xefe be + Xbfb bb ■ (3.39) 

In the absence of nonminimal coupling of CDM (i.e. for universally coupled baryons and 
CDM) the species have the same velocities, i.e. fbbb = febe = f ^ dlndm/dlna. 

In the following we will assume the same baryon-to-CDM ratio for each galaxy and we 
will set this to be the background value, i.e. Xe = Wc and Xb = ujb- However, one could 
also consider different populations of galaxies with different baryon-to-CDM ratios and 
study the effects of equivalence principle violations on large scales between these different 
populations (see e.g. [59]). 
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4 Parametrization 


4.1 Time dependence 

As already mentioned, at the background level the dark energy can be defined by simply 
giving a specific time evolution for the Hubble parameter. For simplicity, we assume that 
the expansion history corresponds to that of rcCDM, so that H is given by 


H\a) = Hl 


n^,oa-^ + (1 - 


(4.1) 


where is the fraction of matter energy density today, tc is a constant parameter and 
the scale factor a is normalized to unity today. This choice of parametrization for the 
background is motivated by the fact that observations suggest that the recent cosmology 
is very close to ACDM, which corresponds to w = —1, and deviations from ACDM in 
the expansion history are usually parametrized in terms of rc 7 ^ — 1. In the absence of 
modihcations of gravity and nonminimal couplings, i.e. for om = 7 c = 0, re coincides 
with the equation of state of dark energy, i.e. u)de in eq. (2.28). Another advantage of 
this parametrization is that the background expansion remains close to the observed one, 
even when ay[ or 7 c are switched on and matter does not scale as a~^ (see eqs. (2.25) 
and (2.26)). In this way we can assume that the background cosmological parameters are 
those fitted by a simple ACDM model. See discussion at the beginning of Sec. 5 and in 
App. A.l. 

In the framework of our effective description, gravitational modifications are encoded 
in the functions a-Q, om and ax; and the non-minimal coupling of CDM is parametrized 
by 7 c-^^ The time dependence of these parameters is undetermined in general. In order 
to obtain some quantitative estimates about how much future observations will be able to 
constrain these parameters, we will focus in the following on a specific functional form for 
their time dependence. 

For simplicity, we will assume that the functions an, om and ax share the same time 
dependence F(t), 


q;b(^) = aB,o r(t), 

aM(i) = om.o r(t), (4.2) 

Q!T(t) = q;x,o r(t), 

where F is normalized to unity today, i.e. F(to) = I, and aB, 0 ) ctM,o and q:x,o denote the 
current values of these parameters, which we wish to constrain. To be more specific, we 

the quasi-static approximation, the parameter uk does not appear in any equation (note that 
the combination c^a does not depend on Qk), while ac,c and ud.c only enter through the combination 
7 c (the combination c^a does not depend on od.c, since Wc = 0), so that their individual values remain 
unconstrained in the analysis. 
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(4.4) 


will consider the following time evolution/^ 




1 - ^m{t) 

1 — rinijO 


where is the total nonrelativistic matter fraction introduced in (2.29) and Om,o hs 
present value. Thus, T vanishes when the unperturbed energy density of dark energy 
is negligible, such as at high redshift, and one recovers general relativity. The above 
parametrization is analogous to the one proposed in [10,14], up to a normalization factor. 

We parametrize the time dependence of 7 c by assuming that the parameter /3.y, defined 
in eq. (3.13), is time-independent, so that 


7c(i) = , 


(4.5) 


and the time dependence on the right-hand side can be computed from eq. (2.14). This 
choice of parametrisation allows to include coupled quintessence [60] as a special case, or 
more generally other cases where the nonminimal coupling of CDM remains active also 
when becomes negligibly small, since one can have = 0 while ^ 0. 

Moreover, vanishes in matter domination, see App. A.2 for details. Therefore, 

when —>■ 1, then T —)• 0 and 7 c —)• 0, which corresponds to the standard matter 

dominated phase for the background evolution. However, while modifications of gravity 
switch off in this limit (i.e. aB)«M)CtT 0), the nonminimal coupling parametrized by 
(3^ remains active (see eq. (4.8) and discussion in the next subsection). 

Let us briefly discuss the theoretical constraints coming from the stability conditions 
[1,3,20]. As discussed in Sec. 2, the absence of ghost-like and gradient instabilities in the 
tensor fluctuations respectively requires > 0—which will be always assumed here and 

in the following—and > 0. Requiring that the second condition is satished at all times, 
eq. (2.5) implies 

q;t,o > “1 • ('^•6) 

For scalar fluctuations, these two conditions become a > 0 and > 0, where the expres¬ 
sions for a and are respectively given in eqs. (2.12) and (2.14). In the following we 
assume that a > 0 is satished by an appropriate choice of the parameters ax, «B and aD,c 
and we will exclude parameters for which the combination c^a (see eq. (A. 5)) becomes 
negative before z = 0 (see again App. A.2 for details). 


4.2 Initial conditions for the perturbations 

We set the initial conditions during matter domination, i.e. when ~ 1, and thus T ~ 0. 
In this limit om — 0 and 7 c ~ 0, so that, according to eqs. (2.22)-(2.23), both CDM 
and baryons behave as conserved species at the background level. Moreover, ckt — 0 and 

Another possible choice would be 

~ o 3™ _L o "I ’ (4-3) 

which has the advantage to be directly related to the scale factor a. We have checked that this choice leads 
to constraints similar to those obtained with the choice (4.4). 
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eqs. (3.12)-(3.13) respectively imply that /3 b — 0 and j3^ ~ 0. Therefore, deep in matter 
domination eqs. (3.15) and (3.16) simplify to 

h + 2HSb 2. [ubSb + ujcSc] , (4.7) 

dc + 2H6c — 2^^ + (l + /3^) Wc<3c] , (4-8) 

where ujb^c are constant. 

This linear system can easily be solved by diagonalizing it. One can find solutions 
written as 

^b ~ ^b,in j — ^c,in <^m) ('^•9) 

with constant and scale-independent bias parameters given by 


bb;m = 


1 -|- /?^Wc — 


2a;2 + (1 - 


2l3^ujc‘^b 


be,in — 


-1 + y^4/32a;2 -h (1 - P^uicY 




7^c 


The respective growth functions Gc and Gb are identical, solutions of the equation 


(4.10) 


G + 2HG - (1 + (3y,be,in) G = 0. (4.11) 

As usual, we will consider only the growing mode solution of this equation, G+. In 
conclusion, we find that baryons and CDM possess spectra that are initially proportional 
and then grow similarly. 

Although we use the full expressions from (4.10) and (4.11) in our numerical analysis, 
it is instructive to consider approximate expressions for small values of /3.y. For small (3-^ 
eq. (4.10) yields 

bb,in = 1 ~ ^“^0^ + C*(/3^) ) be,in = 1 + iOe<-^bl3^ + C){/30 , (4-12) 

while the growing solution of eq. (4.11) is of the form 

G+(a) = + O{(30 . (4.13) 


Thus, for small f3,y the initial conditions in matter domination are simply given by 


Sb{a, A:) ~ (1 - 6 o{k) , 6e{a, fc) ~ (1 uJeUJb00) S^ik) . (4.14) 


4.3 Fiducial models 

For our analysis, we take as fiducial evolution of the Hubble parameter the function 

H{a) = Ho^i}ni,oo>~^ + 1 — I2m,o ) (Fiducial) (4-15) 

which corresponds to the ACDM evolution, i.e. w = —1 in eq. (4.1) and a quantity 
evaluated on the fiducial model is denoted by a hat. The fiducial value for two of the 
parameters that appear in our analysis is taken to be zero, 

om.o = Q;t,o = 0 ) (Fiducial) (4.16) 
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but we consider several options for the parameters (3^ and aB,o- In addition to the simplest 
case where these parameters are zero, it is also instructive to consider fiducial models where 
either of these parameters is nonzero. 

We will distinguish three fiducial models, characterized respectively by the parameters 

I) ACDM: AbjO = P'f — 0; 

II) Braiding: = 0, q;b,o = —0.01, 

III) Interacting: ob.o = 0, = —0.03, 

while the other parameters take the common values prescribed in (4.15) and (4.16). Case 
(I) gives the usual ACDM for the perturbations. In this case the generalized Einstein 
equations and the modified continuity and Euler equations reduce to the standard ones. 
Case (II) corresponds to a mixing between the dark energy and gravity kinetic terms at the 
level of the perturbations. Finally, in case (III) we allow for a non vanishing interaction 
between dark energy and CDM, which is active for perturbations but does not affect the 
background because = 0, and thus 7 c = 0. Let us stress that the background 

evolution is exactly the same for all three fiducial models. 

5 Fisher matrix forecasts 

Our constraints will be based on a Fisher matrix analysis applied to the galaxy and weak 
leasing power spectra [61,62] and to the correlation between the ISW effect in the CMB 
and the galaxy distribution [63]. In general, the Fisher matrix is defined as 

^2 ln£(6») 

de^de^ 

where C is the likelihood function, 0 is a set of parameters. The expectation values are 
over realizations. In the fiducial models I and III 7 c vanishes when varying along j3^ 
(since = 0) and thus, since (3^ = 0 (see eqs. (3.15) and (3.16)), I3.y only appears 

quadratically in the perturbation equations. We have checked that observables depend 
only mildly on 7 c for the fiducial II. Thus, we choose (3"^ rather than [3^ as the independent 
variable in the analysis. In summary, we have the parameters 

6 = {w , aB,o , “M.o 5 aT,o > ■ (5-2) 

Our goal here is to estimate the precision on the above parameters that will be reached 
by forthcoming spectroscopic and photometric redshift surveys with Euclid-like character¬ 
istics [39] (see e.g. [60,64,65] for analogous studies). In particular, we are interested in 
identifying the degeneracies affecting these parameters and their origin. To simplify this 
analysis we will fix the other background cosmological parameters to their Planck esti¬ 
mated values: For w = —1 these are given by [66] h = 0.6731, h‘^^lb,o = 0.0222 and 
h'^Qcfl = 0.1197, while for w ^ —1 we choose the values of Q,b,o Dc,o such as to main¬ 

tain the same angular diameter distance as in the w = —1 case [66]. See details in the 
App. A.l. 
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5.1 Galaxy clustering 

The galaxy power spectrum in redshift space is given by eq. (3.38). Including the correc¬ 
tions due to the Alcock-Paczynski effect, the observed power spectrum reads [67] 

Pohs{z]k,ij) =M{z) [b^{z) + Pmiz,k) , (5.3) 

where the normalization factor M{z) is given by 


N{z) 


H{z)bl{z) 
H{z)D\{z) ’ 


Da{z) 


1 r dz 
^ Z Jq H{z) 


(5.4) 


and Da is the angular diameter distance. Moreover, we assume the bias between galaxies 
and the total matter distribution, 6g = 5^/5ra-, to be scale independent. Its fiducial value 
has little effects on the constraints; in the following we will assume it to be 5g = y/1 + z [68]. 
It can be taken as a nuisance parameter but we will fix it to its hducial value, as a 
consequence of the discussion at the beginning of Sec. 6. Finally, /eff is given in eq. (3.39) 
and Pm{z, k) is the total matter power spectrum, given by 


P^{z,k) = Tliz)Po{k) , 


(5.5) 


where 

Tra{z) = OJbiz) bb,in Gb{z) + Ulc{z) 6c,in Gdz) (5.6) 

is the matter transfer function, Poik) is the initial power spectrum of matter fluctuations, 
(5m, 0 ) during matter domination and bb,in, be,in are defined in eq. (4.10). As the effects 
of dark energy and modihed gravity intervene at late times, the initial spectrum is in¬ 
dependent of the parameters We have neglected corrections due to the shot noise 
in the number of galaxies and the radial smearing due to the redshift uncertainty of the 
spectroscopic galaxy samples and Doppler shift due to the virialized motion of galaxies 
(see e.g. [27,69]), which become relevant on small scales. 

We assume a spectroscopic redshift survey of 15 000 squared degrees, sliced in eight 
equally-populated redshift bins (we take the galaxy distribution as given by [70] with 
a limiting flux placed at 4 x 10“^® ergs“^ cm“^) between z = 0.5 and z = 2.1. The 
corresponding Fisher matrix is given by [62] 


mLSS 

^ab 




V 


bins 


2(27r) 


lirk^dk 


51nPobs(5;; fc,//) 91nPobs(^;^,h) 


'-1 


de° 


deb 


(5.7) 


where V, kmin and femax are, respectively, the comoving volume and the minimum and 
maximum wavenumbers of the bin. In this formula we have neglected the intrinsic statis¬ 
tical error associated with the white shot noise from the Poisson sampling of the density 
field [71]. However, to be conservative, we choose the maximum wavenumber /cmax such 
that the galaxy power spectrum dominates over the shot noise and we are well within 
the linear regime. More specihcally, for each redshift bin we take /cmax as the minimum 

^^Since modifications of gravity affecting the background evolution take place only at late time, we 
are insensitive to the the shift in the matter-radiation equality and to the change in scale of the power 
spectrum turnaround described in [60]. 
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between tt/(2R), where R is chosen such that the r.m.s. linear density fluctuation of the 
matter held in a sphere with radius R is 0.5, and the value of k such that niPg(k) = 1, 
where hi is the number density of galaxies inside the bin. We have checked that these 
values of /cmax are always smaller than + z)), with ag = 400kms“^, i.e. the scale 

where the peculiar velocity of galaxies due to their virialized motion becomes important. 
For the minimum wavenumber, we assume /cmin = Mpc“^. 

Since we work in the quasi-static limit and Po{k) is unaffected by the parameters 0, the 
effects of modifications of gravity and nonminimal couplings are scale-independent. Thus, 
the integration over k in eq. (5.7) simply gives an overall normalisation to the Fisher 
matrix. 


5.2 Weak leasing 


For weak leasing, we consider leasing tomography [72]. The angular cross-correlation 
spectra of the leasing cosmic shear for a set of galaxy redshift distributions ni{z) is given 
by 



Wi{z)Wj{z) 

X^{z) 


ki{z)P^+^[z,ke{z)] , 


(5.8) 


where xiz) = Jq dz/H{z) is the comoving distance and the leasing efficiency in each bin 
is given by 

Wiiz) = xiz)j , (5.9) 

with each galaxy distribution normalized to unity, dzni{z) = 1. Moreover, P$+,i,(A:) is 
the power spectrum of -|- T. Using eq. (3.26), it is related to the matter power spectrum 
by 

P^+'^ik) = T^_^^(z,k)Po{k) , (5.10) 


where 


Tq>+^{z,k) = 


3a^H^ 

2/c2 


[2 -|- CXt + + ^ b ) + /d'y^cbc)\ Tm{z) 


(5.11) 


is the transfer function for T Finally, we define ki{z) = Ijxi.z') as the wavenumber 

which projects into the angular scale t. 

We assume a photometric survey of 15 000 squared degrees in the redshift range 0 < 
z < 2.5, with a redshift uncertainty Oz{z) = 0.05(1-1-2:), and a galaxy distribution [73] 


n{z) oc 2 ^ exp 



(5.12) 


where zq = 2^/1-412 and Zm is the median redshift, assumed to be Zm = 0.9 [27,74]. 
Then, we divide the survey into 8 equally populated redshift bins. For each bin i, we 
define the distribution ni{z) by convolving n{z) with a Gaussian whose dispersion is equal 
to the photometric redshift uncertainty cj^( 2 i), Zi being the center of the ith bin (see 
also [60,65]). 


20 








Neglecting the shot noise error due to the intrinsic ellipticity of galaxies, the Fisher 
matrix for the cross-correlation spectra in eq. (5.8) is given by [75,76] 


pWL 
^ ab 


■^sky 




^Tr 

2 


de^ 



de^ 



(5.13) 


where we choose ^min = 10 and £max = 300. Assuming Euclid-like characteristics [39] for 
the galaxy density and intrinsic ellipticity noise, we have checked that the chosen ^max 
corresponds to scales where the shot noise is negligible and perturbations are only mildly 
beyond the linear regime at small redshift.^^ 


5.3 ISW-Galaxy correlation 

As a third probe, we consider the cross-correlation between the ISW effect of the CMB 
photons and the galaxy distribution in the photometric survey, which is a valuable probe 
of dark energy and of its clustering properties in the late-time universe (see e.g. [77,78]). 
We treat the galaxy survey as for the weak lensing analysis of the previous section, i.e. we 
divide it into 8 bins and, for each bin, we consider the same galaxy distribution. Following 
[79], the projected galaxy overdensity in the bin i is given by 


gi{h)= dz ni{z)bg{z)6in[z,nx{z)] , 


(5.14) 


while the ISW effect is given by 


n = - 


I 


OO Q 

+ ^)[z,nx{z)] . 


(5.15) 


With these definitions, the angular power spectra of the projected galaxy overdensity and 
of the ISW effect are respectively given by 


roo 


ni{z)nj{z)h^^{z) Pra[z,k(,{z)] , 




dz 


/o 


H{z) 

X^{z) 


( dTijy+q, 

V dz 


iz,k)] Poik) 


. k=ki(z) 


(5.16) 

(5.17) 


Analogously, the angular cross-correlation spectrum between the ISW effect and galaxies 
reads 


^ISW-gal^^^ _ _ 


roo \ 

dz { ni{z)bg{z)Traiz) 


10 


X‘^{z) 


dz 


(z,k)Poik) 


- k=kt[z) 


(5.18) 


The Fisher matrix for the ISW-galaxy correlation is given by (see e.g. [80,81]) 


^ ISW-gal _ 
ab 


/sky ^ {2i+l) 


dC 


ISW-gal 


f=f 

C — -Cji 


de° 




ISW-gal 


(^) 


de^ 


(5.19) 


^^Notice that the value of £max chosen here is smaller than what is usually assumed in comparable 
analyses (see e.g. [27] and references therein). 
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where we use £min = 10 and £max = 300 and the covariance matrix is given by 

Covjk{i) = , ( 5 . 20 ) 

where is the full CMB angular power spectrum. We have omitted from this 

expression the CMB noise, which is negligible for CMB experiments such as WMAP and 
Planck, and the galaxy shot noise. We have checked that the latter is small up to the 
chosen £max- 

6 Results 

In this section we present the results of the Fisher matrix analysis and the associated 
degeneracies between parameters. We start by discussing the effects of nonstandard gravity 
on the evolution of homogeneous quantities. As shown below, they are important to 
understand the effects on perturbations. 

6.1 Background 
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Figure 1: Relative change of the baryon and CDM density fractions, with respect to their fiducial 
values, as a function of the redshift z, depending on the values of the parameters re, ob.O) cpvi.O) 
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Before presenting the results of the Fisher matrix analysis, we discuss how the back¬ 
ground evolution is modihed when one goes slightly away from any of the fidneial models 
by modifying one of the parameters. The results are summarized in Fig. 1, where we have 
plotted the evolution of the difference between and their respective hducial value. 

As is clear from (2.25), the parameter fib is only affected by a change of the background 
history embodied by H{z) or by a variation of the effective Planck mass M. It is thus only 
sensitive to a change of the parameters w or om- In the former case, the evolution of pb, 
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Fid. 

Obs. 

10^ X (t(1 -I- w) 

10^ X cr(aB.o) 

10^ X cr(aM,o) 

10^ X (T(aT,o) 

104 X a(/32) 

I 

GC 

7.0 

18.6 

24.5 

- 

1.4 

WL 

1.6 

4.3 

42.1 

- 

5.7 

ISW-g 

15.5 

4.4 

20.2 

- 

31.3 

Comb 

1.6 

3.0 

14.6 

- 

1.35 

II 

GC 

7.2 

18.6 

33.8 

24.4 

2.7 

WL 

1.4 

4.4 

67.4 

98.9 

6.4 

ISW-g 

5.0 

4.2 

24.5 

43.2 

56.0 

Comb 

1.3 

3.0 

19.0 

20.8 

2.5 

III 

GC 

0.22 

0.40 

0.22 

0.22 

1.4 

WL 

0.17 

2.12 

0.18 

0.18 

5.7 

ISW-g 

0.88 

2.78 

0.88 

0.87 

31.3 

Comb 

0.13 

0.39 

0.14 

0.14 

1.4 


Table 1: 68% confidence level (CL) errors on each individual parameter, assuming that the others 
take their fiducial values, for each fiducial model and observable: galaxy clustering (GC), weak 
leasing (WL), ISW-galaxy correlation (ISW-g) and the combination of the three (Comb).^^ The 
parameter aT,o is unconstrained in fiducial model I, see explanation in Sec. 6.2.1. 


and thus Life, is modified because H is changed. In the latter case, the evolution of pb does 
not change but that of does. These changes are independent of the other parameters 
and one does not need to distinguish between the three fiducial models. 

For VLc, the situation is exactly the same as Q-b when w or ccm are changed, provided 
there is no coupling between dark energy and CDM, i.e. 7c = 0. This is apparent in the 
boxes corresponding to the fiducial models I and II, for which (3^ = 0. By contrast, if 
we start from the fiducial model III, where (3^ ^ Q, and modify either w or q;m; then 
the deviation of He with respect to its fiducial value is amplified due to the coupling 
7 c generated by a nonzero CgO^I'^ combined with a nonzero /3.y. For the same reason, 
i.e. 7c / 0, we observe a deviation of He when q;t,o or ap.o are switched on, in contrast 
with the other fiducial models. This also explains why one sees a deviation from the 
fiducial model II when (3^ is switched on. 

The modifications of the background quantities discussed above affect the observables 
both indirectly, through their effect on the evolution of perturbations, and directly, because 
the observables explicitly depend on H and Hm (see for instance eq. (5.11)). Therefore, a 
qualitative analysis of the effect of the parameters 6 on the observables is rather complex 
and must take into account both the background evolution and the quantities T^^c and 
Tiens- This is why we resort to a Fisher matrix analysis, which allows us to quantify the 
combined effects on the observables. 
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Figure 2: Modifications of the evolution of perturbations from their fiducial values, as a function 
of redshift, for the different parameters w, ae.O) cim.O) cit,o and f]~f. From top to bottom, relative 
variation of the effective growth factor /es, eq. (3.39), the matter transfer function T^, eq. (5.6), 
the Weyl potential transfer function T, 5 +ii<, eq. (5.11) and its derivative with respect to redshift, 
for the three different fiducial models (respectively 1, II and III, from left to right). As 
^zT^+^^r vanishes in matter domination, we have normalized it to its value at z = 0 instead of its 
value as a function of the redshift. 


6.2 Forecasts 

Let us now discuss the results of the Fisher matrix analysis. The unmarginalized errors on 
the parameters are summarized in Tab. 1 while the two-dimensional contours are presented 
in Figs. 3, 4 and 5. Red dotted, green dashed and yellow solid lines respectively correspond 
to galaxy clustering, weak lensing and ISW-galaxy observables. The combination of the 
three observables, given by summing the three Fisher matrices, is plotted in thick solid 
black line. The shaded blue regions in the plots correspond to instability regions, where 

^'^Our constraints on are in qualitative agreement with those obtained for coupled quintessence 
in [60], taking into account that the parameter defined in this reference is related to ours by = 2/3^. 
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c^a < 0 .^^ 

For each observable, the Fisher matrix including all the parameters is ill-conditioned 
and cannot be inverted. This means that the observables do not have the constraining 
power to resolve the degeneracies (see e.g. [82]). Thus, when plotting the two-dimensional 
contours we do not marginalise over the other parameters but we hx them to their fiducial 
values. 

As shown in Tab. 1, the forecasted constraints from the three probes for the same 
fiducial model are comparable, within an order of magnitude. This reflects the comparable 
effects on the observables, shown in Fig. 2, given our choice of A^max and £max for the 
spectroscopic and photometric surveys, respectively, which translates into a comparable 
number of modes for the three probes. More precisely, the effects of gravity modihcations 
and nonminimal couplings is slightly larger on the lensing potential and ISW effect but 
this is compensated by a larger number of modes in the spectroscopic survey. 

Specihcally, for this survey the number of modes is roughly given by Amodes ~ Abins x 
V X (47r/3)(/cmax/27r)^, where Abms = 8 is the number of bins and V is the (average) 
comoving volume of the bins. Assuming A:max = 0.1/iMpc“^, this yields Amodes ~ 10®. 
For the photometric survey we have Abodes ~ Abins x /s^y x ~ 3 x 10®. As a rule 
of thumb, the relative effects of aB ,05 om,o and ot.o on the three observables are typically 
of the order of 0(0.1) at redshift z ~ 1, see Fig. 2. Thus, one expects to be able to 
constrain these parameters at the level of 0(0.1)“^ x i.e. few percents (which is 

improved by one order of magnitude for fiducial III, where the effects on the observables 
are larger), if all the other parameters are hxed. The ISW-galaxy correlation is limited 
by cosmic variance but due to the larger sensitivity of to the modifications of 

gravity, it sometimes provides constraints comparable to those from the other probes.^® 
The effect of /3^ is typically of the order of a few at redshift z ~ 1 and this parameter can 
be constrained at a level of a few x 10“^ for galaxy clustering and weak lensing. Given the 
smaller effect on the ISW and the smaller number of modes for the photometric survey, 
the ISW-galaxy correlation provides the weakest constraints on this parameter. We also 
notice that the degeneracy of this parameter with the others is rather small. 

6.2.1 Fiducial I: ACDM 

Let us study the constraining power of the observables around a ACDM model. The errors 
are reported in Tab. 1 and the 68% CL contours are shown in Fig. 3. In Tab. 2 we report, 
for each Fisher matrix, the eigenvector associated to the maximal eigenvalue (called here 
maximal eigenvector), which provides the direction maximally constrained in parameter 
space, i.e. the one that minimizes the degeneracy between parameters. 

At first view, the parameter ot.o seems to contribute to the growth of perturbations 
through the combinations T/, and Tc, defined in (3.18), and to the lensing potential through 
the combination Tiensj given in (3.27). However, it turns out that these combinations in 
fact do not depend on ut for this choice of hducial model. 

^®Here we conservatively exclude the instability region from the allowed parameter space. A more 
refined treatment would require multiplying the likelihood function by a theoretical prior that excludes 
the forbidden region, which is impossible to achieve with a Fisher matrix analysis (our priors cannot be 
represented with an invertible matrix). 

^®We thank Alessandro Manzotti and Scott Dodelson for pointing out a numerical underestimation of 
the noise in the ISW-galaxy correlation in an earlier version of this paper, corrected here. 
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Obs. 

Fiducial I 

Fiducial II 

Fiducial III 

GC 

(0.012,-0.007,0.005,0,1) 

(0.022,-0.013,0.007,0.01,1) 

(-0.626,0.348, -0.629,0.64,1) 

WL 

(-0.345,-0.115,-0.007,0,1) 

(-0.463, -0.136, -0.001,0.004,1) 

(1, -0.074,0.910, -0.914, -0.293) 

ISW-g 

(0.053,0.7,0.154,0,1) 

(0.856,1,0.117,0.063,-0.609) 

(-0.997, -0.138, -0.989,1, -0.068) 

Comb. 

(-0.008,-0.012,0.005,0,1) 

(-0.055, -0.034,0.006,0.009,1) 

(1, -0.285,0.953, -0.964, -0.867) 


Table 2: First eigenvector of the Fisher matrices, for the basis {w, aB, 0 j c«m, 0 i Q^t.Oj with 
the maximum eigenvalue, corresponding to the combinations of parameters that are maximally 
constrained by experiments. The coefficients are normalized by the maximum component and 
rounded to three significant digits. 


More precisely, when w = —1 and = 0, one finds that 


'^b,c = ax + /3| = ctt + CgO = —2(1 + aB)C + — 3^ . (6.1) 

^ da H 


When one goes away from the fiducial model by switching on the parameter ap, while all 
the other parameters keep their fiducial value, one gets /3| = —ot so that the dependence 
on ax vanishes in It is immediate to check that ax disappears in Tiens for the same 
reason. Thus, the parameter ax,o cannot be constrained by a Fisher matrix analysis for this 
choice of fiducial and will be dropped from the analysis in this subsection. Correspondingly, 
the component in the ax,o direction of the maximal eigenvectors vanishes, see Tab. 2. 

Let us now examine the situation when ae is switched on while all the other parameters 
take their fiducial value. The T combinations are then given by 


Tb,c = 




Tiens 


c^a ’ 


( 6 . 2 ) 


with 

c^a =—(2 + 3nin)aB — 2a|. (6.3) 

For small values of a-Q, we thus find 


T 


6,C 


2 + 3flr 


-aB , 


Tiens ^ - 


2 + 3Qr 


-Ob • 


(6.4) 


Thus, one expects the impact of as to increase as flm diminishes, which is in agreement 
with the results plotted in Fig. 2. 

When one changes om from its fiducial value (the other parameters keeping their 
fiducial value), one finds 

Tft^c = Tiens = O'M ■ (6.5) 

As seen in Fig. 2, the effect of au and a-Q on the growth of structures (i.e. on /eir and 
Tjn) is roughly the same in magnitude but opposite in sign, which is in agreement with 
the relations found in (6.4) and (6.5). This qualitatively explains the degeneracy observed 
in the aB,o“CkM,o panel of Fig. 3 for galaxy clustering and the corresponding components 
of the maximal eigenvectors in Tab. 2. By contrast, the degeneracy between ob and om 
observed for weak leasing does not seem to agree with the values of Tiens in (6-4) and (6.5). 
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FIDUCIAL I 



lO^x (1+w) l0^xaB,o lO^x (I'M,0 


Figure 3: Two-dimensional 68% CL contours for the fiducial model I (ACDM model), obtained 
by fixing all the other parameters to their fiducial values. The parameter ot.o is absent, as it is 
unconstrained on this fiducial model. Shaded blue regions correspond to theoretically forbidden 
parameter space where c^a < 0. Note that the axis range is different for different parameter planes. 


The reason for this discrepancy is that the background is also modified when ayi ^ 0, as 
discussed earlier, whereas the background for ap 0 is the same as the fiducial one. 
Since the transfer function depends not only on the coefficient Tiens but also on the 

background, the degeneracy is more complex. In fact, the background modification also 
affects the matter growth but more modestly than for weak lensing. 

To conclude, let us note that a large region of the observationally constrained parameter 
space is forbidden by the stability requirements, i.e. c^a > 0. 
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FIDUCIAL II 



10^x(l+w) lO^x Q-flo lO^x QfA /,0 lO^x aro 


Figure 4: Two-dimensional 68% CL contours for the fiducial model II (braiding model with 03,0 = 
—O.OI), obtained by fixing all the other parameters to their fiducial values. Shaded blue regions 
correspond to c^a < 0. The axis range is different for different parameter planes. 


6.2.2 Fiducial II: Braiding 

For this fiducial model, we have the value q;b,o = —0.01, where the negative sign is to 
satisfy the stability conditions. This corresponds to dark energy models where the kinetic 
term of vr comes from a mixing with gravity [4, 6] , which are sometimes called braiding 
models [83,84]. The unmarginalized errors are reported in Tab. 1 and the 68% CL contours 
are shown in Fig. 4. Note that the allowed parameter space is much larger than in the 
previous fiducial because for 03,0 7^ 0 the null energy condition can be violated without 
instabilities [4]. 
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In this case, c and Tjens depend on ckt: their partial derivatives with respect to ot 
on the fiducial model are given by 

^ 9nl QTiens ^ - 2 - 2aB) . , 

daT ( 3 flm + 2 + 2 aB)^ ’ daT (Sflm + 2 + 203)^ ’ 

which confirms that this parameter must be included in the analysis. 

For this fiducial, the plane q:b^o~«t,o in Fig- 4 has the same background evolution as 
ACDM. Therefore, all the effects are controlled by and Tiensj so that the degeneracies 
can in principle be understood analytically from their expressions in terms of 03,0 and 
aT,o- For instance, for small 03,0 and ax.o one finds 


T 


b,c 


3aB,o — 1) (2aB,o + (2 — 3nm) aT,o) 

as.o (6flm + 4) + 4aT,o 


{1-nm) (0.54aT,o - 0.6AaB,o) , (6.7) 


where in the last equality we have expanded at linear order for small 1 — flm and used 
q;b,o = —0.01 + Aq!b,o- This explains the degeneracy between Aq;b,o and q;t,o observed 
in the growth. By the same procedure we find Tiens — (1 — fim) (O.lSax.o — l-2AaB^o)) 
which explains why Aob.o is more constrained than q;t,o by lensing observations. 

Similarly to fiducial I, the effect of changing 03,0 and aM,o on the growth of structures 
is roughly the same in magnitude and opposite in sign. This effect can be qualitatively 
understood by expanding T;, ^ for small Aan^o and aM, 0 ) analogously to what was done in 
Sec. 6.2.1. This degeneracy cannot be seen for the lensing, because the modifications of 
the background also play a role. 


6.2.3 Fiducial III: Interacting 

In this model we have a nonzero fiducial value for the parameter (3^ = —0.03), which 

implies an active coupling between CDM and dark energy. The unmarginalized errors 
are reported in Tab. 1 and the 68% CL contours are shown in Fig. 5. Notice that the 
constraints for this fiducial model are generally stronger than those for models I and II 
(see below). As one can verify in Fig. 2, this is due to the enhancement of the effects on 
the observables, caused by the nonminimal coupling. 

In this case, q;t,o must be included in the analysis, because and Tiens depend on 
q;t,o through the term (3^(3^. Indeed, let us examine the case when aT,o and um.o are 
switched on while w = —1 and 03^0 = 0. Using CsO?-!'^ = y^2(aM — cit) = —\/2/3^ (we 
assume um > ckt to satisfy the stability condition) one finds 

Tfe = UM - \/aM - ttT , Tc = UM - -v/aM - OT + Wc^c) + P'if (6.8) 

and 

Tiens = — \/Um “ «T (^■ybcUJc ■ (6.9) 

However, the degeneracies observed in Fig. 5, for example in the plane aM, 0 “«T, 0 ) cannot 
be understood directly from the above expressions because, as we saw in Fig. 1, the 
background is modified, not only when ayi (or w) is changed but also when ax is changed. 

Another notable degeneracy appearing in Fig. 5 is between w and the parameters 
—Q!x,o or aM,o- This can be partially understood from the fact that w appears in the 
combination 

c^a ~ 3(1 + w){l - Hm) - 2(aM - ot) = 3(1 - Um) (1 + rc - um.o + aT,o) , (6.10) 
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FIDUCIAL III 



Figure 5: Two-dimensional 68% CL contours for the fiducial model III (interacting model with 
/3-y = —0.03), obtained by fixing all the other parameters to their hducial values. Shaded blue 
regions correspond to c^a < 0. The axis range is different for different parameter planes. 


where we have used r/ ~ —u)(l — Llm) in eq. (A.4) for the first equality and rim.o — 1/3 in 
the last one. However, background effects play an important role as well. 

The term in eqs. (3.18) and (3.27) translates here as —— ctT id-y, see eq. (6.8). 
This term encodes the new effects that arise when both modifications of gravity and 
nonminimal couplings are considered, as emphasized in [20]. These effects explain the 
qualitative difference, in the size and shape, between the contours of fiducial III (Fig. 5) 
and those of the other two fiducial models. Not only are the constraints tighter by an 
order of magnitude in this case, but also the maximal eigenvectors of the Fisher matrices 
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point in different directions, see Tab. 2. 


7 Summary and conclusions 

In this paper, we have investigated the consequences of both modifying gravity and allow¬ 
ing a coupling between CDM and dark energy. If the propagation speed of dark energy 
is not too small, one can rely on the quasi-static approximation because the small scale 
fluctuations of dark energy have the time to relax to the quasi-static regime [38]. In 
this case, the parameters describing deviations from ACDM, which are usually four for 
Horndeski-like theories [3,9,10], reduce to three: a-Q-, am and ax [20]. Moreover, the 
coupling of a fluid of CDM particles conformally and disformally coupled to dark energy, 
can be described by a single parameter 7 c, see eq. (2.24). 

The dynamics of matter perturbations also simplifies. In particular, as discussed in 
Sec. 3, it is described by a system of two coupled equations, eqs. (3.15) and (3.16), re¬ 
spectively for baryons and CDM. In these equations, the four parameters above enter in 
three combinations (see eq. (3.18)): ax, (3^ (a combination of ob, om and ax) and (3^, 
the latter describing the nonminimal coupling of CDM perturbations. As explained in 
more details in Sec. 3, these distinctions are frame-dependent, as one can verify using the 
relations (2.10) (see also [20] for more details). 

The growth of fluctuations is usually described in terms of the growth rate, which 
modulates the galaxy power spectrum in redshift space and can thus be measured with 
redshift space distortions. We have computed the effective growth rate for galaxies made 
of baryons and nonminimally coupled CDM, in the presence of modifications of gravity. 
This is the first general treatment of this kind, to our knowledge. 

Deviations from the ACDM model can also affect the propagation of light through their 
effect on the scalar Weyl potential, i.e. the sum of the two metric potentials in Newtonian 
gauge. A fourth parameter, (3b (proportional to ob), together with the three parameters 
above, is necessary to fully describe this effect, which can be measured in the weak leasing 
and ISW effect (see eq. (3.27)). 

As discussed in Sec. 4, the evolution of perturbations depends on the time dependence 
of the Hubble rate and of the parameters described above. In the present work we have 
taken the Hubble rate to be the same as in tcCDM. Moreover, the parameters aB, om 
and ax grow as 1 — Dm, so that modifications of gravity disappear in matter domination, 
while the nonminimal coupling remains active at all times, i.e. (3^ = constant. We have 
studied the constraining power of a future redshift survey with Euclid specifications on 
the parameters w, aB,o, ckm.o and ax,o- 

More specifically, in Sec. 5 we computed the Fisher matrix of the galaxy power spec¬ 
trum, the weak leasing power spectrum as well as the correlation spectrum between the 
ISW effect and the galaxy distribution. We have considered five parameters, namely w 
(describing the background evolution), the current values of aB, aM and ax, and the con¬ 
stant nonminimal coupling parameter /3^, and assumed three fiducial models: (I) ACDM, 
(H) a braiding model with aB,o = —0.01 and (HI) an interacting model with f3.y = —0.03. 

The unmarginalized 68% CL errors on these parameters are reported in Tab. 1 in 
Sec. 6. For the current values of aB, an and ax, the errors are of the order of 10“^-10“^ 
for fiducial models I and H and an order of magnitude better for the fiducial model HI. The 
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error on /3^ is of the order of 10“^ for all fiducial models. Given the large number of free 
parameters and the degeneracies among them, the Fisher matrices cannot be inverted to 
compute the marginalized contours. Therefore, we have shown the two-dimensional 68% 
CL contours in Figs. 3, 4 and 5 —together with the excluded parameter space from stability 
conditions—by setting all the other free parameters to their fiducial values. Moreover, we 
have provided a discussion on the origin of the degeneracies and the constrained directions 
in parameter space in Tab. 2. As shown by the contour plots, all the three observational 
probes are complementary in breaking degeneracies in parameter space. 

This analysis can be generalized in several directions. First, the background cosmolog¬ 
ical parameters should be included in the analysis as nuisance parameters. In this case, it 
is important to take as well into account other cosmological data such as the CMB, the 
baryon acoustic oscillations and the supernovae Type la. Another direction is exploring 
alternative parametrizations of the background evolution and/or of the time dependence 
of the parameters ob, oM) cut and (3^. For instance, assuming that the a’s vanish at 
early times, as we did, considerably limits the effect of dark energy on certain observables 
such as the CMB or the matter power spectrum. On the other hand, one could assume 
other equally motivated time dependencies (even different for different parameters), which 
are expected to lead to larger effects in the observables. The final goal is to extend this 
analysis beyond the quasi-static approximation to include larger scales and other species, 
such as neutrinos and photons. Such a program has been initiated with the development 
of the publicly available Boltzmann codes EFTCAMB [85] (see [86] for a recent applica¬ 
tion to Horava gravity) and COOP [87]. In this case, at least one more parameter, ok, 
must be considered in the analysis. On the other hand, one may expect that some of the 
degeneracies found in this paper can be resolved by the scale dependence appearing once 
the full dark energy dynamics is taken into account. We leave this for future work. 

Acknowledgements. We wish to thank Tessa Baker, Phil Bull, Scott Dodelson, Pedro 
Ferreira, Zhiqi Huang, Alessandro Manzotti, Valeria Pettorino, Federico Piazza, Ignacy 
Sawicki and Emiliano Sefusatti for stimulating discussions. J.G., M.M. and F.V. thank the 
APC {AstroParticule et Cosmologie) and PCCP (Paris Center for Cosmological Physics) 
for kind hospitality. F.V. acknowledges partial support from the grant ANR-12-BS05-0002 
of the French Agence Nationale de la Recherche. 

A Details on the parametrization 

In this appendix we provide some details about the determination of the background 
parameters in our numerical calculations and about the value of the effective functions in 
our parametrization. 

A.l Background quantities 

Assuming that gravity is standard at recombination, dark energy can only affect the 
best fit value of the cosmological parameters inferred through the measurement of the 
comoving distance to last scattering with the CMB spectrum. Thus, we assume that the 
comoving distance to last scattering is fixed and given by its best fit measurement [66] 
and we compute the values of the background cosmological parameters inferred from this 
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observation. Let us discuss how these are determined. When w = —1, these are chosen as 
the base ACDM best fit values of the Planck TT+lowP parameters [66]. When w ^ —1, 
we determine the initial conditions for the background matter components by requiring 
the comoving distance^' 

r 

XiZin;klrafi,w) = dzH~^ {z] klrafl, w) (A.l) 

Jo 

to be the same as the one of the ACDM model. More precisely, for each value of w, we 
associate the parameter Qrn,o{w) defined by the relation 

XiZin] klm,o{w),w) = xizin] = “1) (A.2) 

where we have on the right hand side the standard ACDM value, evaluated by using the 
value ^ = 0.02222/i-2 + 0.1197/i-2, with h = 0.6731, which 

corresponds to the estimate deduced from the measurements by the Planck satellite [66]. 
We take Zin = 100, deep in the matter dominated era, when the effects of dark energy are 
negligible. 


A.2 The combination 


Here we provide details on the calculation of For convenience we define the pa¬ 

rameter 


^ - 3 r ^ ^ + (1 - L!„,o)a-3- 


(A.3) 


which enters naturally in eqs. (2.25) and (2.26). For om = 7 c = 0, the fraction that 
appears on the right hand side reduces to the energy density fraction of dark energy, 
1 — Dm, but this is not the case in general. From eq. (2.14), the combination c^a reads 


cla = (1 + aB)(3 -3ri- 2() - 3Dm - 2^ , (A.4) 

where rj and ^ are defined above, respectively in eqs. (A.3) and (3.12). By using eqs. (4.2) 
and (4.4) and the background evolution equations (2.25) and (2.26) to evaluate Ab in this 
expression, this can be written as 


c^a = 3(1 - Dm - 7)+aB 


l-3rj{l + 2 




1-D, 


- 2(aM - 37c Wc) 


Dc 


1 - Dr 


— 2Q:g — 2 q!t(1 4“ Ob) + + ob) • 

Finally, one can replace Xc by its expression (4.5) given in terms of 

The equation (A.5) is thus a quadratic equation for X = CsO^^'^ ^ of the form 


(A.5) 


- HA - C = 0, 


(A.6) 


where 

B = (A.7) 

i i ijn 

^^The comoving distance is related to the luminosity distance and the angular-diameter distance 
Da by the relations Dl{z) = (1 + z)x(z) and Da{z) = x( 2 )/(l + z). 
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and 


C = 3(l-0m 


r?) + 


1 - 3?7 


1 + Oin 
1 


2 

1-^^m 


Om 


aB~2aB —2aT(l+Q!B)^+2Q;M(l+«B) • 

(A. 8 ) 


Let us extract from this quadratic equation the relevant solution. 
Let us start with the case q:b,o = 0, which implies 


B = 0, C = 3(1 - Om - ?7) - 2 (q;t - om) (aB,o = 0) 
in (A. 6 ), and the solution is therefore 

, (aB,o = 0). 


(A.9) 


(A.IO) 


Both signs of this solution can be chosen and lead to the same phenomenology as long as 
the sign of 13^ is chosen to obtain the same 7 c. In the matter dominated era, corresponding 
to rim —>■ 1 and 7 —)• 0 , the stability condition thus imposes 


aM,o > «T,o («B,o = 0) • 


(A.ll) 


Let us now consider the case aB,o / 0. In the past limit flm —^ 1, one finds that 


B = V2(3, 


'^OJc 


rimOB 

l-rim 


V2(3, 


'^LOc 


aB,o 


1 - n 


m,0 


(A.12) 


behaves like a constant, while C —)• 0. Consequently, the two solutions of the quadratic 
equation in this limit are A = 0 and X = \/2/3.^a;caB,o/(l ~ rim,o)- la order to recover a 
standard matter dominated regime with 7 —)• 0 in the past limit flm Ij one needs to 
pick up the A = 0 solution in the past. This determines the choice of the sign among the 
two solutions 


A = 


B±^^WTW 


(A.13) 


which yield A = (B ± |B|)/2 in the limit flm —^ 1- One thus concludes that, depending 
on the sign of aBfl P'y, the solution is 




Cga 


1/2 


B - + 4C 

- j -, 

±VC , «B,0 

b + Vb^T^ 

2 ■ 


iP'y OB.O > 0) 

= 0 ) 

i(3'y a-Qfl < 0) 


(A.14) 

(A.15) 

(A.16) 


As above, both signs on the right hand side of eq. (A.15) can be chosen. The stability 
condition c^a > 0 is obtained by requiring that the above solutions are real. 


B Matter evolution equations in a generic frame 

For completeness, we provide here the evolution equations for matter in a generic frame. 
In a generic frame where both baryons and CDM are nonminimally coupled, eqs. (3.15) 
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and (3.16) read 


+ (2 + 3jt)H6b = + Th)6m , (B.l) 

+ (2 + 3^c)H6c = + Te)<5m , (B.2) 

with 

Tf, = ax + /?! + + 2/3.yj,/3^)a;b6fc + + 13^(13^^ + 13^J] iojjc , (B.3) 

Tc = ax + /3| + [f3-ibf3^c + l3^iP'yb + /^7c)] ^bh + (/3^^ + 2(3^^(3^)ojchc ■ (B.4) 

For the case discussed in the main text of minimally coupled baryons, i.e. j3^^ = 0, one 
recovers the expressions in eq. (3.18). 


Under a frame transformation (2.8), a;/ = does not change. Moreover, in the 

quasi-static limit the density contrasts 5i does not change either (the explicit transforma¬ 
tions are discussed in [20]). In particular, this implies that bj = bj. Therefore, by using 
the transformations of the a’s given in eq. (2.10), 77 given in eq. (3.19) and those of the 
/I’s given in eq. (3.20), one finds the expressions for in the frame g^u, 

Tb = (1-I-Tft)(l-I-ao) - 1, 

(B.5) 

Te = (l + Te)(l + aD)-l. 

Using the expressions above and that the factors 1 -|- od is cancelled by the change of time 
between the two frames, di = sJC/{I + aD)dt, one can check that the form of eqs. (B.l) 
and (B.2) is frame-independent. 
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